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ABSTRACT 

Context. The amount of molecular gas is a key for understanding the future star formation in a galaxy. Because H2 is difficult to 
observe directly in dense and cold clouds, tracers like the CO molecule are used. However, at low metallicities especially, CO only 
traces the shielded interiors of the clouds. In this context, millimeter dust emission can be used as a tracer to unveil the total dense 
gas masses. However, the comparison of masses deduced from the continuum SIMBA 1.2 mm emission and virial masses (thought to 
trace the whole potential of the clouds) in a sample of giant molecular clouds in the Small Magellanic Cloud (SMC) by Rubio et al. 
(2004) and Bot et al. (2007) showed a discrepancy between these two quantities that is in need of an explanation. 
Aims. This study aims at better assessing possible uncertainties on the dust emission observed in the sample of giant molecular clouds 
from the SMC. It focuses on the mass comparison in the densest parts of the giant molecular clouds where CO is detected in order to 
confirm the mass discrepancy observed. 

Methods. New observations of the south-west region of the SMC were obtained with the LABOCA camera on the APEX telescope. 
All the giant molecular clouds previously observed in CO are detected and their emission at 870/mi is compared to ancillary data. The 
different contributions to the sub-millimeter emission are estimated, as well as dust properties (temperatures, emissivities), in order to 
deduce molecular cloud masses precisely. 

Results. The (sub-)millimeter emission observed in the giant molecular clouds in the south west region of the SMC is dominated by 
dust emission and masses are deduced for the part of each cloud where CO is detected and compared to the virial masses. The mass 
discrepancy between both methods is confirmed at 870/jm with the LABOCA observations: the virial masses are on average 4 times 
smaller than the masses of dense gas deduced from dust emission, contrary to what is observed for equivalent clouds in our Galaxy. 
Conclusions. At present, the origin of this mass discrepancy in the SMC remains unkown. The direct interpretation of this effect is 
that the CO linewidth used to compute virial masses do not measure the full velocity distribution of the gas. Geometrical effects and 
uncertainties on the dust properties are also discussed. 

Key words. ISM: clouds — submillimeter — ISM: molecules — Magellanic Clouds 



1. Introduction 

Star formation is observed to occur within molecular clouds. The 
star formation rate in a galaxy therefore depends on the amount 
of molecular gas available, but this quantity remains difficult to 
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estimate precisely. In particular, at low metallicities, the rela- 
tionship between H2 and CO (the most widely used tracer of 
molecular gas) remains unclear. 

The Small Magellanic Cloud (SMC) is one of the closest 
and easiest to observe, low metallicity galaxies. Due to its prox- 
imity, molecular clouds in this galaxy can be resolved at vari- 
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ous wavelength. Numerous molec ular clouds in the Magel lanic 
clouds w ere ob served in CO lines dlsrael et al.ll9 93: Rubio ^tafl 
ll993alH fl996) and were found to have weak CO emission. If the 
CO molecule is tracing accurately the molecular gas, then the 
star formation per unit CO emission in the SMC is higher than 
in most galaxies. Suc h an effect is obse rved in other low metal- 
licity dwarf galaxies dLerov et al.l l2006) and could be an "evolu- 
tionary effect" or the fact that CO observations largely underes- 
timate the amount of molecular gas in these conditions. Indeed, 
at low metallicities, CO can be photodissociated and trace only 
the densest parts of the clouds, while H2 molecules have the ca- 
pacity to self-shield and would be present at larger radius than 
the CO molecule. 

This "underestimation" by CO of the molecular content 
of SMC molecular clouds has been seen and studied mainly 
through the Xco = Nh 2 /Ico factor. Using CO data only, the 
velocity dispersion observed by the CO line width is thought 
to trace the gravitational potential of the cloud and therefore 
the total molecular content, which can be compared to the 
CO intensity in the same region. Analy zing NANTEN data 
for the whole SMC, iMizuno etafl d2001l) fou nd X co = (1 



5) x 10 21 cm _2 (K.km.s" 1 )- 1 . Blit z et alJ d2007l) corrected these 
data for sensitivity and resolution effects and found Xco - 
(1 - 1.5) x 10 21 cm _2 (K.km.s -1 ) -1 . These values are systemat- 
ically larger than the canonical galactic value Xco - 2.3 x 
10 20 cm _2 (K.km.s _I )~ 1 dStrong et al.ll 19881) . supporting the idea 
that CO is present only in part of the molecular cloud. There 
is evidence for a d ependence of the Xco factor with radius 
dRubio et alJ Il993ah and at the smallest resolved scales Xco 
determined for individu al resolved clouds might become close 
to the galactic va lue dBolatto et al.1 l2003t llsrael et al.l 12003b 
iBolatto et alll2008l) . 

To better understand the total molecular content and its rela- 
tionship to CO emission, other tracer of the dense and molec- 
ular gas at low metallicities have been looked for. Dust con- 
tinuum emission can be used to trace the dense and cold inter- 
stellar medium where the gas is molecular and it can the refore 
be us ed to unveil the total mass of dense gas in clouds. Ilsraell 
d 19971) was the first to introduce this method in the Magellanic 
Clouds, comparing IRAS and CO data, and he deduced higher 
H2 masses than wh at was found using CO alone. More recently, 
iLerov et all (l2007h attempted to spatially trace molecular hydro- 
gen using far infrared emission observed with Spitzer data in the 
SMC. In dense regions, they observed large quantities of dust not 
associated to either Hi or CO emission. A study of t he star form- 
ing re gion N83 confirms this result on smaller scale dLerov et al.l 
2009). However, the dust emission in the far-infrared is sensi- 
tive to the dust temperature, leading to large uncertainties in the 
cloud masses. In the sub-millimeter and the millimeter range, 
dust emission is difficult to observe but de pends only l i nearly 
on the dust temperature and is optically thin. lRubio et al.l d2004) 
observed for the first time with SIMBA on the SEST telescope, 
the 1 .2 mm emission from dust in a quiescent molecular cloud 
in the SMC. The dense gas mass deduced were several times 
higher than th e dynamical ma ss deduced from CO observations 
in this cloud. iBot et all d2007) extended this study with SIMBA 
to a larger sample of molecular clouds in the south-west region 
of the SMC. In all clouds, the gas mass deduced from the dust 
emission is systematically larger than the virial mass deduced 
from CO, even for conservative values of the dust emissivities 
and free-free contribution to the millimeter emission. All these 
studies come to agreement in the fact that molecular clouds in 
the SMC are more massive than what can be deduced from CO 



emission. This mass difference is understood in a scenario where 
CO traces only the densest parts of the molecular clouds. 

However, in case of virial equilibrium, it has been as- 
sumed that the observed motions should balance the gravita- 
tional pressure and the dynamical masses obtained from molec- 
ular line observations should trace the total molecular masse s 
(e.g.lMacLaren et al.|[l988l:lRosolowskv et akEool iDu & Yan 
20081: iDunham et al.ll20ia) . The result from IBot et al.l (120071 
where the virial masses in the SMC systematically underesti- 
mate the true mass of the clouds (as traced by the dust mass) 
remains unexpected. Partial support of the clouds by a mag- 
netic field has been invoked to explain this mass discrepancy. 
However, these results need to be confirmed with observations in 
the sub-millimeter range in order to more certainly exclude any 
possible contamination to the fluxes (free-free emission, CO line 
contribution), as well as unknown instrumental effects and dust 
properties. This study focuses on the parts of the giant molecular 
clouds that are detected in CO. These regions represent only the 
highest density parts of the giant molecular clouds but are most 
appropriate for the comparison with CO and the check of a mass 
discrepancy. 

We present new observations of the sample of molecular 
clouds in the south west region of the SMC taken with the 
LABOCA camera on the APEX telescope (section|2|i. These ob- 
servations at 870/mi complement the SIMBA data at 1.2mm, 
CO and radio data, as well as 160/mi Sp itzer data from the 
combined S 3 MC an d SAGE-SMC surveys dBolatto et al.ll2007b 
iGordon et al.l 12009!) . We build spectral energy distributions of 
each giant molecular cloud (GMC) detected both in the sub- 
millimeter and in CO and assess the possible contaminations and 
uncertainties(Sec. l2.3b . The dust emission in each cloud is then 
used to determine dense gas masses for each molecular cloud 
in the south-west region of the SMC. These mass estimates are 
compared to virial masses and confirm the mass discrepancy 
(Sec. [3). 



2. Data and methods 

2.1. The data 

The dust continuum observations were obtained using the L arge 
APEX Bolometer Camera (LABOCA, ISiringo et all 120091) on 
the Atacama Pathfinder Experiment ( APEX. iGtisten et alJ 2006) 
telescope in Chile. LABOCA is a bolometer array consisting of 
295 receivers packed in a hexagonal structure resulting in a to- 
tal field of view of 1 1 'A. It has a central frequency of 345 GHz 
(~ 870yum) with a band width of 60 GHz and the beam-size 
(HPBW) was measured to be 1972. 

The 870/mi LABOCA observations were done during 
October 2007. The weather conditions were mostly good to ex- 
cellent with a precipitable water vapor (PWV) content typically 
between 0.5 to 0.9 mm resulting in zenith opacities of the or- 
der of 0.2. The pointing was checked regularly on strong contin- 
uum sources and planets, and the pointing accuracy was found 
to be mostly better than 3". The focus settings were determined 
from observations of planets typically once per night and during 
changes of the weather conditions, sunset and sunrise. To obtain 
the atmospheric opacity the sky emission was measured regu- 
larly about every two hours. Calibration factors were derived by 
observing Neptune and Uranus as well as secondary calibrators. 
Relative offsets and relative gains (flatfield) of the individual 
bolometer channels are determined about once per month and 
are provided by the APEX staff. The scientific mapping was car- 
ried out in spiral raster mode with the raster pattern alternating 
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between 2x2 and 3x3 spirals to achieve an optimal sampling of 
the region of interest. Nine positions raster maps were done in a 
140" spiral at l"/s speed and 60° inclination. Raster lengths were 
varied in x and y directions between 1000" and 1200", and the 
spacings between 500" and 600". Dithering of 100" to 120" from 
the center position was performed. Additional observations of 
the central 6' x 6' region in a four by four points raster were per- 
formed to increase the signal to noise ratio. In total we obtained 
101 maps and a pure on-source integration time of 8.9 hours. 
The final rms achieved was 8mJy/beam. 

The LABOCA data were reduced using the reduction pack- 
age BoAQ (Schuller et al. in prep). The processing of the raw 
time series, i.e. the time-ordered data stream (down-sampled to 
25 Hz) of each channel and scan, consists of the following steps: 
calibration correction by applying a linearly interpolated cali- 
bration factor; opacity correction by applying the linearly in- 
terpolated opacity at the elevation of the scan; correction for 
temperature drifts due to the cryosystem using two bolometers 
that have been sealed to block the sky signal for this purpose; 
flat fielding by applying the relative bolometer gains; conversion 
from counts to Jansky; flagging of bad channels. In addition to 
the known dead or noisy channels, bad channels were identi- 
fied as those having an rms of their time series a given ratio r 
higher or 1/r lower than the median rms of all channels. The fur- 
ther reduction applied (in part on several occasions) the follow- 
ing tasks: correlated noise removal on the full array; correlated 
noise removal on groups sharing the same amplifier box and the 
same wiring; baseline subtraction; despiking; flagging of data 
outside suitable telescope scanning velocity and/or acceleration 
limits; flagging of bad channels. Each reduced scan (i.e. raster 
pattern) was then gridded into a weighted intensity map and a 
corresponding weight map. The weights of the data points con- 
tributing to a certain pixel of the intensity map were determined 
as 1 1 a 1 where cr denotes the rms of the reduced time series of 
the corresponding channel and subscan. Individual maps were 
then coadded, again noise-weighted, to build the final intensity 
map and the corresponding weight map, which in turn allows to 
retrieve the rms for each pixel and to construct a signal-to-noise 
(S/N) map. 

The described data reduction is affected by the presence of 
astronomical signals in the time series, which leads to a flux loss 
due the subtraction of correlated noise. In case of strong sources 
this is manifestly revealed by the appearance of areas of negative 
fluxes adjacent to the sources. Moreover, the source emission 
wrongly contributes to the determination of the rms and weights. 
To minimise these effects we applied an iterative approach by 
subsequently improving a model of the flux distribution of the 
astronomical source. 

In a first step a map was produced by a "blind" execution of 
the reduction as described above. A first source model was con- 
structed from the resulting map by extracting all pixels above a 
given S/N threshold (typically 3) and setting the remaining pix- 
els to zero. This model was converted to a time series for each 
channel and subtracted from the databefore any baseline sub- 
traction, despiking and correlated noise suppression. After the 
weight determination and before the map was built, the model 
was re-introduced into the time stream. The resulting map was 
used to extract a new source model. In subsequent iterations 
this process is repeated until the measured flux distribution con- 
verges. As a result, most of the faint extended emission should 
be recovered in the reduction. 
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The 20' x20' image of the south-west region of the SMC that 
was obtained, is presented in Fig.Q] It is the first sub-millimeter 
map of a large region in the Small Magellanic Cloud. This south- 
west region of the SMC has had most of its known molecu- 
lar clouds mapped at a 43" resolutio n and at hight sen sitivity 
in the CO (J=l-0) line with SEST (iRubio et al J 1 1993atfbh . All 
the GMCs observed in CO lines with SEST (shown as con- 
tours in Fig. [TJ are detected with the present LABOCA obser- 
vation. There is also extended emission well outside observed 
CO peaks. If the 870^/m emission originates from dust, this re- 
sult shows that there is much more dense gas than what is traced 
by current CO observations. Comparing with the MIPS 160/im 
map of the same region, we see that the spatial distribution of 
the 870/im emission is well correlated with the one observed in 
the far-infrared. A short analysis of the dust emission outside the 
observed CO is presented in Appendix iBl 

The present paper compares the 870/im and CO data and is 
therefore restricted to the regions where CO has been observed 
and detected. The entities for which we deduce masses in the 
following and that are refered to as "clouds" are the densest parts 
of the molecular clouds, associated with CO emission. 

For comparison, we complement these data with the S IMBA 
observations at 1 .2 mm previously studied dBot et all l2007h 
and with Spitzer MIPS observations at 160//m in the same 
region from the combined S 3 MC and SAGE-SMC surveys 
(iBolatto et all 120071; iGordon et at] 120091) . Furthermore, in or- 
der to estimate the contribution of free-free emission to 
the observed millimeter emission, we also use radio ATCA 
maps at 8.6, 4.8 GHz that were obtained from J. Dickel 
(http://www.phys.unm.edu/~johnd/). 



2.2. Flux estimates 

To estimate the LABOCA, SIMBA, MIPS and radio fluxes in a 
way that is comparable to CO, we do the following. All maps 
are convolved to the 43" CO(1-0) resolution and projected on 
the same sampling grid. Because we are interested in compar- 
ing the LABOCA, SIMBA, MIPS and CO fluxes, it is important 
to remove the extended emission in all maps in order to com- 
pare emission in the region detected with all tracers (the peaks 
of the GMCs). To do so, we begin by masking the region where 
CO is detected and fit in each map (LABOCA, radio, SIMBA, 
MIPS) a 2D plane to the emission outside the CO detected re- 
gion. The fitted plane is then removed from the maps (in units 
of surface brightness) to leave only the peak of emission where 
CO was observed. The integrated fluxes are then computed in 
a region above a defined Wco threshold. This threshold is de- 
fined individually for each cloud so that the area of the cloud we 
sample is similar to the one that was used to determine the viria l 
masses and CO luminosities of the clouds ( Rubi o et alJ 1993b). 
This method enables us to analyze the emission coming from 
the dust associated with the CO emission only. It attempts to fil- 
ter extended envelopes around the clouds - uncertainties on this 
filtering process will be discussed in section l4TI - as well as pos- 
sible emission asso ciated with Hi. Using an H i map of the Small 
Magellanic Cloud dStanimirovic et aLlH999l) . we estimated that 
the remaining Hi column densities associated to the peaks of the 
molecular clouds we study, after filtering, are less than a few per- 
cent of the total hydrogen column densities as determined from 
the dust millimeter emission. The fluxes obtained are summa- 
rized in Tab. [T]and are used to build the spectral energy distribu- 
tion for each cloud (c.f Fig. [2j. 
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Table 1. Obse rved fluxes from the far-i nfrared to the radio for our sample of giant molecular clouds. The radii and velocity disper- 
sions AV from Rubio et al.l(ll993bl) (and lLequeux et al] (I 19941) for Hodge 15) used to compute the virial masses are also mentioned. 
Notes: (1) Association of the molecular clouds [RLB93] SMC-B2 1, [RLB93] SMC-B2 2 and [RLB93] SMC-B2 6; (2) Association 
of the molecular clouds [RLB93] SMC-B2 3 and [RLB93] SMC-B2 4. 



name 


SIMBAD name 


Radius 


A.V 


S 160//m 


S 870//m 


S l. 2mm 


S 3cm 


S 6cm 






pc 


km.s~' 


Jy 


mJy 


mJy 




mJy 


mJy 


LIRS49 


LI-SMC 49 


16.7 


6.2 


16.1 ± 1.5 


1358 ± 113 


403 


±66 


28 ±3 


6±3 


LIRS36 


LI-SMC 36 


18.6 


3.8 


7.2 + 0.8 


1015 ± 104 


258 


±45 


16 ±4 


20 ±3 


SMCB1#1 


[RLB93] SMC-B1 1 


13.8 


3.1 


1.6 ±0.4 


234 ±4 


83 ± 


25 


< 6 


< 5 


SMCB1#2 


[RLB93] SMC-B1 


16.2 


4.9 


8.9 ± 2.0 


789 ± 146 


175 


±76 


15 ±3 


14 ±2 


SMCB1#3 


[RLB93] SMC-B1 3 


12.9 


3.2 


2.9+1.4 


265 ± 84 


59 ± 


19 


12 ±5 


13 ±3 


Hodge 15 


[H74] 15 


18 


4.6 


2.1 ±0.8 


751 ± 176 


245 


±89 


< 9 


10 ±3 


SMCB2 South 1 




15.4, 10.9, 15.0 


3.7,4.7,4.5 


18.8 ±8.9 


1660 ± 500 


422 


±211 


114 ±23 


37 ± 10 


SMCB2 North 2 




19.8, 11.4 


4.6, 2.9 


21.5 ± 5.1 


1505 ± 280 


278 


± 125 


40 ± 10 


24 ±6 


SMCB2#5 


[RLB93] SMC-B2 3 


15.0 


4.2 


< 3.7 


172 ± 116 


57 ± 


44 


< 9 


4±2 



2.3. Origin of the 870fim emission 

At long wavelength like 870/im, the broad-band flux densities 
determined from the maps may contain non-negligible contri- 
butions by thermal free-free continuum emission and CO line 
emission. 

2.3.1. Free-free contribution 

In iBot et al] (I2007I) . estimates of the free-free emission and of 
the CO line at 1.2mm were found to be negligible contributions 
to the measured fluxes. Even if we expect this to be the case at 
870^m, it is safe to quantify the effect here too. Furthermore, 
the free-free estimates obtained in IBot etalJ (120071) were derived 
from the Payn e et alj (120041) radio catalog, which had integrated 
fluxes in regions that could be different in size with respect to 
the molecular clouds we are sampling, and when no source was 
present in the catalog it was not clear whether this was a non- 
detection or a non-observation. In this section, the use of the 
radio maps directly enables us to derive much more accurate es- 
timates of the free-free contribution to the (sub-)millimeter emis- 
sion for the regions we are interested in. 

All regions observed at 870//m and in CO(1-0) also show 
radio emission that is detected at 4.8 or 8.64 GHz. A few regions 
(e.g. LIRS49) show an inverted spectrum (which could be due 
to the different spatial filtering between the observations at the 
two frequencies). But in most cases, the radio spectral index is 
flat, consistent with an origin of the emission that is mostly free- 
free. Indeed, for free-free emiss ion, the flux den sity scales as 
0.95 - 0.16 In v with v in GHz (lRevnoldslll992l) . and is fitted 
to the radio emission detected at the highest frequency (plain 
black line in Fig. 0. Comparing this scaling with the emission 
at 870pm and 1.2 mm, the free-free contribution to the (sub- 
millimeter emission is not significant (we estimate it to be 1 
to 8% of the 870pm emission and remove it from the observed 
LABOCA and SIMBA fluxes). 

2.3.2. CO(3-2) line contribution 

The J=3-2 CO emission contribution to the 870pm broad band 
emission is co mputed using the m easured CO(3-2) line, when 
available, from Nikolic et al. <l2007h . or otherwise, by applying a 
high CO(3-2)/CO(2- 1 ) ratic0 of 2 to the CO(2-l) observed emis- 
sion. Using a Jy/K conversion factor of 8.4 and a bandwidth 
for the LABOCA bandpass of 60GHz, we find that the CO(3- 



2 we take the ratio of the velocity integrated intensities in K.km.s 



2) line contribution to the observed emission is insignificant in 
all clouds (it represents at most 0.3% of the emission at 870pm). 

The 870pm and 1 .2mm fluxes in each cloud are then domi- 
nated by dust continuum emission. 

2.4. Method 

In Sec. 12.3. Il and l2.3.2l we showed that the emission detected by 
LABOCA at 870pm is dominated by dust emission. The method 
to determine hydrogen masses fro m dus t emission is then similar 
to the one described bv lBot et"aL1 (120071) . In short, the hydrogen 
column density in each molecular cloud can be deduced from: 



e 87() fi87()(7ciHjr) 

where Tdust is the dust temperature, e^ () is the emissivity of dust 
per hydrogen atom at 870pm and B \ is the Planck law. In order to 
deduce dense gas masses from the dust emission detected with 
LABOCA, we therefore have to know the dust temperature in 
the clouds of the SMC and the dust emissivity at 870pm. 

2.4.1. Dust temperatures and emissivity index 

In Bot et al.l (l2007l) . it was assumed that the spectral dependance 
of this dust emission in the (sub-)millimeter is similar to the one 
in our galaxy (S v oc v^B v (r £ /„ s/ ) with /3 = 2.0) with a single 
temperature of 15 K, typical of GMCs in the Milky Way. Even 
though the dust temperature dependance of the submillimeter 
emission is mild, the uncertainty on this parameter can create 
scatter in the mass determined for the sample. Both the dust 
temperature and the emissivity index in the molecular clouds of 
the SMC can be better constrained by comparing the LABOCA, 
SIMBA and Spitzer MIPS fluxes for each cloud. 

A modified blackbody fits to the MIPS, LABOCA and 
SIMBA fluxes is done for each region with both the dust tem- 
perature and the emissivity index as free parameters (dashed 
line in Fig 0. Such a fit gives values of ft between 1.9 and 
4 (c.f. Tab. 13 and very low dust temperatures. However, high 
values (e.g. ft ~ 3) are not believed to be physical and could 
be a consequence of filtering in the SIMBA data (c.f. 0, even 
though our focus on emission peaks attempts to limit this effect. 
Alternatively, it could be due to a mix of temperatures within the 
regions we sample. 

It should be noted that since there is most probably a range 
of temperatures within the beam, the dust temperature deter- 
mined from a single modified black body fit is an upper limit 
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Fig. 1. Top image: LABOCA image in the south west region of 
the SMC (Units are in Jy/beam). Middle and bottom images: 
LABOCA images convolved to a ~ 40" resolution with CO and 
MIPS 160/mi contours overlaid respectively. All the molecular 
clouds observed in CO are detected at 870pm. Extended emis- 
sion is well recovered in the LABOCA map and the 870pm 
emission matches well the 160/mi one. 



Table 2. Results from modified black body fits to the MIPS, 
LABOCA and SIMBA fluxes for each molecular cloud. The 
two first values correspond to a fit where both the emissivity 
index and the dust temeprature are free to vary, while the two 
last columsn corresponds to the dust temperatures determined 
for fixed values of the dust emissivity index (J3 = 2 and /3 — 1). 



name 


P 




dust 


T I3=2 
dust 






dust 




LIRS49 


2.5 = 


t0.2 


11.2 + 0.3 


13.1 


+ 


0.2 


20.6 


±0.4 


LIRS36 


3.3 = 


t0.2 


8.5+0.3 


12.0 


+ 


0.2 


18.0 


±0.4 


SMCB1#1 


1.9 = 


t0.3 


12.2 ±0.8 


11.7 


+ 


0.3 


17.1 


±0.7 


SMCB1#2 


3.8 = 


b0.4 


8.3 + 0.5 


13.1 


+ 


0.4 


20.9 


± 1.0 


SMCB1#3 


3.8 = 


t0.5 


8.2 + 0.6 


13.5 


+ 


0.9 


22.6 


±2.2 


Hodge 15 


2.4 = 


t0.5 


9.2 + 0.7 


10.3 


+ 


0.4 


14.1 


±0.7 


SMCB2 South 


3.2 = 


b0.6 


9.3 + 0.9 


13.1 


+ 


0.8 


21 ± 


2 


SMCB2 North 


4.0 = 


t0.3 


7.8 + 0.3 


13.1 


+ 


0.4 


21.0 


±0.9 


SMCB2#5 


2.2 = 


t 1.0 


11.7 + 2.4 


12.1 


+ 


1.3 


18 ± 


3 



to the mean temperature. Since lower temperatures would give 
higher dust masses, our simple fit gives lower limits on the mass 
and is thus a conservative estimate to determine the dense cloud 
masses. 

Using the spectral energy distributions built for each cloud 
(Fig. [2]) we observe that in all the molecular clouds of this study, 
the emission is consistent with a modified black body with a 
standard spectral index of — 2. Fixing the dust emissivity to 
this standard value, we observe low dust temperatures (12 K on 
average) for all the molecular clouds studied here. The temper- 
atures range is small (from 10.2 to 13K) and consistent with a 
constant temperature at a 3cr level. However, these low dust tem- 
peratures (which correspond to a radiation field intensity that 
would be 10 times lower than the cano nical value from the so- 
lar neighbourhood (Mathis et al. 1983)) are surprising. Indeed, 
some regions like LIRS36 or LIRS49 harbour active star form- 
ing regions and our regions could contain, at least partly, dust 
heated by the surrounding stars. However, because our study fo- 
cuses on regions where CO is detected, it is possible that the 
dust in these regions is predominantly cold. Even though, it is 
surprising that dust temperature in giant molecular clouds of the 
SMC are lower than what is observed for molecular clouds in the 
Milky Way, especially since at low metallicity the dust shielding 
will be reduced. 

The temperature and emissivity index of the dust associated 
with the 870pm emission is therefore still difficult to constrain 
and more data in the 100-300pm wavelength regime would be 
needed to do so. In the following, we assume j3 — 2 and the 
temperatures (T^ust ~ 12K) determined in such case from the 
fitting procedure, but we also compute lower limits on the mass 
estimates using a low value of the emissivity index: /3 — 1 . A fit 
to the spectral energy distributions with f3 = 1 (dashed blue line 
in Fig. [2]) leads to dust temperatures of Tdust ~ (last column 
in Tab.|2]). In the following, we will therefore also compute lower 
limits on the molecular cloud masses using Tdust — 19 K. 

2.4.2. Dust emissivity at 870pm 

To determine gas masses from the dust (sub-)millimeter emis- 
sion at 870pm it is crucial to know the dust emissivity in this 
wavelength regime. This quantity can be expressed as follow: 

e^ = K &70 x d ium H (2) 

where k^o is the absorption coefficient per unit of dust mass, Xd 
is the dust to gas mass ratio (in the following, we adopt a dust to 
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Fig. 2. Spectral energy distribution from the far-infrared to the radio, observed for each cloud. The fluxes observed at 160, 870/im 
and 1 .2mm are fitted with a single modified black body either with an emissivity index of 1 (dashed blue line) or 2 (plain red line) 
or with a free emissivity index (dotted green line). Radio fluxes (or upper limits) computed in the same regions are overplotted and 
the extrapolated free-free emission is displayed as a black line. This shows clearly that the sub-millimeter and millimeter emission 
is dominated by dust emission rather than free-free. 



gas ratio 0. 1 7 times the one in the solar neighbourhood, as de- 
termined bv lBot et~a"D ((2007) from abundance studies) and [imu 
is the gas weight per hydrogen, taking into account the contribu- 
tion of He. 

The long wavelength emissivity of dust in the dif u se me dium 
has been modeled and tabulated by iLi & Draind (1200 ll) . At 
850^m, they give /csso^m = 0.4cm 2 /g. However, this value may 
not be representative of the dust properties in molecular clouds. 
In particular, coagulation processes can occur between grains. 



The large fluffy aggregates created that way have different ab- 
sorption/emission properties than those in the diffuse medium. 
Enhancements of the dust emissivity can be attributed to coagu- 
lation proc esses. Such an emissivity enhancement is observed in 
our g alaxy dCambresv et alJl200ltlStepnik et aIll2003tlBot et all 
2007) and for giant molecular clouds, the dust emissivity in the 
(sub-)millimeter regime is enhanced by a factor 2-3 with respect 
to the diffuse medium. 
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In order to determine the dust opacity appropriate for the 
GMCs probed by our observations, we chose to estimate it di- 
rectly from sub-millimeter observations in a Galactic molecular 
environment. Without further knowledge, we make the assump- 
tion that the molecular ring in our Galaxy can be used as a ref- 
erence for the giant molecular clouds of the SMC. In particular, 
this assumes that dust evolution (like coagulation) occur in the 
same way in SMC GMCs as in our galaxy. 

The method used is similar to the one bv lBot et"ai] (120071) : a 
map of the dust emission in our Galaxy at 870//m is created using 
FIRAS data. In the molecular ring, the 870/iin dust emission cor- 
related with Hi is substracted from the observed emission using 
the correlation observed at high galactic l atitude. The rema ining 
emission correlates with the CO intensity (Dam e et al . 2001) and 
a linear fit on the correlation gives: 



^870, 



W C0 



'"" - 0.251 +0.003MJy/sr.(K.km/s)- 1 . 



(3) 



We checked that in individual GMCs in our solar neighbour- 
hood, this ratio is similar (~ 0.26 + 0.18- the uncertainty here 
reflects the dispersion between clouds) and is therefore char- 
acteristic of the dust emission in GMCs in the Milky Way. 
Assuming a standard conversion factor between the CO intensity 
and the molecular gas column densities (Xco _ N(H2)/Wco _ 
1.8- 1 20 mol.cnr 2 (K. km/s) 1 dDame et alj200U iGrenier et all 
20051)). a dust temperature in galactic GMCs of 15 K as in 
Bot et al. (2007), we can deduce a dust emissivity at 870/itn for 
the molecular ring. By multiplying this value with the SMC dust 
to gas ratio, we obtain a dust emissivity suitable for SMC GMCs: 



4iQimi H 2' SMC ) = ( 3 - 94 ± °- 05 ) X KT 27 ar 



1 .cm 2 . 



(4) 



This emissivity of dust in molecular regions corresponds to an 
opacity k&jo = 1.26 + 0.02cm 2 /g. 

The uncertainties quoted above for the different values cor- 
respond to formal uncertainties from the fitting procedures and 
do not reflect the intrinsic scatter in the 870//m-Wco correlation, 
nor the systematic uncertainty due to the assumptions which are 
by far the dominant source of error. For example, it is unclear 
whether "coagulation" of dust grains occurs in the SMC and this 
will be discussed further in Sec. 14.21 However, within our cur- 
rent knowledge, most systematic effects will increase the mass 
estimated from the dust emission. 

Using eq.|4]and a dust temperature between 12 and 19 K, we 
can therefore determine gas masses directly from the observed 
870yum dust emission continuum in the CO detected GMCs of 
the SMC. 



3. Results 

Gas masses for all the GMCs of the south west region of the 
SMC observed in CO are deduced from the LABOCA sub- 
millimeter emission and are summarized in Tab. [3] These masses 
are compared to the masses deduced from the application of the 
virial theorem on the CO line width and rad ii (c.f. Tab [3]and Fig. 
|3). For comparison purposes, we took from iBotet all (120071) the 
masses determined from dust millimeter emission and the virial 
masses for a reference sample of similar molecular clouds in the 
solar neighbourhood. 

For self-gravitating clouds, the velocity dispersion of the CO 
clumps balances the gravitational pressure and the virial masses 
should be similar to the cloud masses deduced from dust emis- 
sion. However, we observe that the masses of the GMCs as 
deduced from dust emission are systematically larger than the 



Table 3. mass estimates for the individual GMCs observed in 
the south-west region of the SMC. The first two columns give the 
mass estimates computed from the LABOCA 870/mi emission 
for dust temperatures of 12 K and 19 K respectively, while the 
third column lists the virial masses obtained from the CO line 
widths and the cloud radii. 



name 


M™"'"(12K) 
10 4 M o 


M™" m (19 K) 
10 4 M o 


M yir 


LIRS49 


56 ±5 


31 


+ 3 


12.2 


+ 0.5 


LIRS36 


50 ±5 


23 


+ 2 


5.1 


±0.5 


SMCB1#1 


12 ±2 


5 i 


: 1 


2.5 


±0.2 


SMCB1#2 


33 ±7 


18 


±4 


7.4 


± 1.4 


SMCB1#3 


11 ±4 


6d 


: 2 


2.5 


±0.2 


Hodge 15 


50+12 


17 


±4 


7.2 


±0.4 


SMCB2 South 


65 + 22 


35 h 


: 12 


14.4 


±2.3 


SMCB2 North 


85 ± 11 


46 


±6 


9.8 


± 1.9 


SMCB2#5 


8±6 


4 d 


:3 


5.0 


±0.3 



10.0 



1.0 



0.1 



Galactic clouds (Bot et al. 2007) x 



SMC clouds 
SMC clouds 



this study, T dl 



= 12K • 



this study, T dust =19K 



Ho0ge15 
L!rs36 I 



10 4 



10 3 



Fig. 3. Comparison between the two molecular mass estimates: 
masses deduced from sub-millimeter dust emission and virial 
masses. The mass deduced from the LABOCA emission in the 
SMC are shown for a dust temperature as estimated from the 
SED fitting with /3 = 2 (Td us , ~ 12K, black diamonds) and for 
a dust temperature of 19 K as an upper limit (green rhombs). 
Galactic GMCs studied with the same method are shown for 
comparison (blue stars). In the SMC, the molecular mass de- 
duced from the sub-millimeter are systematically larger than the 
one deduced from the virial theorem applied on CO data, even 
with an upper limit on the dust temperature, in opposition to 
what is seen for equivalent galactic clouds. 



virial masses determined for the same clouds. Our study with 
LABOCA data at 870/i/m therefore confirms the results obtained 
bv lRubioetaLl (12001 and lBot et all (120071) with SIMBA data. 
We find a median mass ratio M vir /M™ m - 0.21 (0.41 for the 
upper limit on the dust temperature), contrasting with the mass 
ratio obtained for equivalent clouds in the Milky Way, which is 
always above 1 . 
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The most straightforward interpretation of the mass discrep- 
ancy we observe in the SMC is then that the CO linewidth do not 
measure the full velocity distribution of the gas. There are two 
independent reasons for why this could be true. 

The CO emitting regions may be present only in the interior 
of the mo lecular clouds, in regions shielded from the UV radia- 
tion field dLequeux et al.lll994tlGio"ver & Mac Lowl2010h . If the 
(unresolved) CO emitting clumps are not evenly distributed in 
the region detected by the CO observations, the observed line 
width might not reflect the full velocity dispersion in the cloud. 

The CO line could come fr om the densest regio ns created by 
shocks in a turbulent medium dKlessen et alj|2009l) . Simulations 
in a turbulent interstellar medium show that these dense re- 
gions with CO mi ght not always be the most shielded ones 
dGlover et al. I l2o!oh . In this case, the velocity dispersion in the 
CO emitting regions would not be representative of the velocity 
dispersion of the molecular gas located also in lower density re- 
gions. The CO line width would then underestimate the velocity 
dispersion of the molecular gas due to its preferential location 
in the densest and less turbulent regions. Such a difference has 
been observed between the CO and HI line width in high lat- 
itude clouds in the Solar Neighbourhood on a few parsec scale 
(Barr iault et al.l20Toh . but it is unclear if such a difference would 
exist on larger scales where many shocks at different locations 
of the cloud must be occuring. 

Due to the low metallicity of the SMC, a Galactic analog to 
understand the CO emitting clouds in the SMC may be denser 
molecular environments. Indeed, if the H2 abundance varies as 
the product of the mean gas number density and the metallic- 
ity dGlover & Mac Lowl feOlO). the motions observed in CO and 
the mass observed in the dust continuum can be compared to the 
velocity of molecular cores and the relatio n to their envelo ppes 
in our Galaxy. In "extinction super cores". iKirk et alJ (120071) ob- 
served that the inner N2H + or C ls O starless cores show a velocity 
dispersion that is systematically lower than the one expected for 
a gravitational support of the region in extinction. The motions 
between the densest molecular entities are hence not sufficient 
to support the enveloppes against gravity, similarly to what we 
observe in the giant molecular clouds in the south-west of the 
SMC. 



4. Discussion 

Our analysis of the two mass estimates was done in a way to 
avoid as much as possible the known caveats (CO-devoid en- 
velope, enhanced dust emissivity): fluxes are extracted from re- 
gions restricted to where CO is observed, we attempt to remove 
the extended emission surrounding the CO peaks and the dust 
emissivity is taken for molecular environments. The impact of 
such effects on the masses deduced is checked in Appendix lAl 
Here we discuss the possible remaining biasses that could ham- 
per the analysis and result in the mass discrepancy. We also dis- 
cuss alternative interpretations for the mass discrepancy we ob- 
serve. 

4. 1 . Unresolved CO clouds and geometrical effects 

In the scenario where CO is confined inside a large H2 envelope, 
then the dust emission will trace the total cloud (CO-emitting 
cloud+envelope). In order to minimize the effect of the CO- 
devoid envelope on the dust emission, we performed a local 
estimate of the extended emission around the CO-detected re- 
gion and removed it from the measured (sub-)millimeter fluxes. 



However, the effectiveness of this removal depends on the size 
of the envelope and on the density profile of the cloud. 

For the same configuration, the virial mass is obtained from 
the CO velocities and the radius of the CO cloud. This radius 
may be overestimated from the observations if the CO clumps 
are not resolved and the virial masses will be overestimated. The 
dimension differences between the dust emitting region, the CO 
emitting region and the area deduced from the observations (i.e. 
the aperture) may then still bias the ratio M,,„/M mm in either di- 
rection (enhancing or decreasing the mass ratio). 

The exact magnitude of such biases is a complicated func- 
tion of the real structure of the cloud (e.g. density), the defini- 
tion of the CO cloud and aperture, and the extent of any CO-free 
dusty envelope surrounding the cloud. We tested the magnitude 
of a geometric bias for a range of density profiles with a simple 
model. In this model, the CO cloud and the envelope are spheri- 
cal. A circular aperture defines the radius that enters the virial 
mass calculation, while also defining the impact of the back- 
ground subtraction to the mass deduced from the dust emission. 
We find that for a simple model in which dust emerges from a 
cloud twice the size of the CO cloud and the aperture covers 
the portion of the cloud not showing CO emission, the bias al- 
most never exceeds ~3Q%. However, if the CO cloud is actually 
substantially smaller than the values quoted in Table [T](e.g. CO 
clumps are not resolved at the SEST resolution), if the dust en- 
velope is very large compared to the CO cloud, or if the aperture 
extends too far from the surface of the cloud, these effects may 
be larger. Conversely, for shallow density profiles and aperture 
definitions too close to the surface of the cloud, it is possible to 
over-correct for a dusty envelope, meaning that our measurement 
of M mm would in fact be too low and exacerbating the difference 
between M VJ > and M mm . This effect has been explor ed in more 
detail s for the N83 region in the Wing of the SMC dLerov et al.l 
2009) and could play a role in the mass discrepancy despite our 
attempt to limit the mass comparison to the CO peaks of the 
molecular clouds. 

To further assess these biases, we performed the mass calcu- 
lations without the extended emission subtraction and with de- 
convolved CO radii. The results are presented in AppendixlAl It 
shows that, as far as can be tested it with the current observa- 
tions, both effects seem to enhance the mass discrepancy. 

Observations with better resolution and comparisons to mod- 
els and simulations of cloud structure should offer the chance to 
better constrain whether geometric effects can drive the observed 
difference. Here we note the possible concern and that we have 
done our best to correct for it given the data available (via the 
extended emission fitting procedure). 

4.2. The dust to gas ratio and the dust emissivity 

In order to compute molecular cloud masses from the dust 
emission, we assumed that the dust properties in the Small 
Magellanic Cloud are similar to the one we observe in our 
Galaxy, except for the lower dust content with respect to the gas 
due to the metallicity. Hence, one possibility to understand the 
mass discrepancy we observe is to question our understanding 
of dust emission in the SMC. 

Unfortunately, our current knowledge of dust grains in the 
SMC is limited and the dust emiss ion in the SMC cou ld be dif- 
ferent than in our Galaxy. Indeed, llsrael et al] (|2010|) observed 
a millimeter emission excess in the Small Magellanic Cloud, on 
the scale of the galaxy, t hat is not explained by current standard 
models of dust emission dBot et al1l2010l) . The nature of this ex- 
cess remains unclear but might be associated to spinning dust 
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emission. Given the fact that dust emission in the SMC is not 
well understood on a global scale, the interpretation of a simi- 
lar excess emission on the scale of individual molecular clouds 
is therefore precarious. We acknowledge this uncertainty but do 
our best within the current knowledge on dust properties in the 
SMC, i.e. we assume that large dust grains at thermal equilib- 
rium do most of the sub-millimeter dust emission. 

Biases due to our uncertainties on the dust temperature and 
emissivity index (note that these quantities are not independent) 
were assessed through the use of an extreme value j3 — 1 giving 
Tdust ~ 19 K for a fit of the observed far-infrared to millimeter 
SEDs. In Fig. [3] we observe that the lower limit on the molecular 
cloud masses obtained with these parameters are still larger than 
the virial masses. 

In Sec. 12.4.21 we assumed that the dust emissivity is equiva- 
lent in the SMC to what is observed in GMCs of the Milky Way, 
and scaled it accordingly to the dust-to-gas ratio difference be- 
tween the two galaxies. However, there is no ground based evi- 
dence that this is the case. First, the knowledge of the dust-to-gas 
ratio in the SMC is quite uncertain. Second, the grain opacity 
could be different. In the following, we discuss how these two 
quantities could vary. 

The dust-to-gas ratio estimate w e use is based on abundance 
measures toward the star AV 304 (Rolles ton et alj 20031). us - 
ing galactic depletions computed from Sofia & Mever (200 1blal) . 
These abundances are consist ent with all HII region abun- 
dances published for the SM C ( Vermeii & van der Hulsti r2002; 
iPerez-Montero & Dfazl2 005: Lebouteiller et al. 2008) expect for 
one (N13, for which a close to solar O/H abundance has been 
deduced) so that we believe our metallicity estimate is accurate 
for the SMC. The depletion factors could be higher than in the 
Milky Way, enhancing the dust-to-gas ratio and decreasing the 
mass discrepancy. However, in order to compensate the mass 
discrepancy we observe with a change of the dust-to-gas ratio, 
all metals would have to be depleted onto dust grains. It there- 
fore seems unlikely that a change of the dust-to-gas ratio could 
explain the mass discrepancy we observe. 

Grain coagulatio n can make the sub -millimeter dust emissiv- 
ity larger than the iLi & Draind (UpOl) value referenced for the 
diffuse medium. An enhancement (a factor of 3) is observed in 
Galactic GMCs. We assumed that the same enhancement factor 
applies to the SMC clouds in our study. Could this enhancement 
be larger in the SMC and explain the mass discrepancy? Taking 
an upper limit on the dust temperature, we observed in the last 
section that the virial mass for the SMC GMCs is at least 2 to 
5 times smaller than the mass deduced from the dust millime- 
ter emission. In order to get M v! > = M"™, the dust emissivity 
in the SMC would then need to be at least 2 to 5 times larger 

than the value we use. We note that such an enhancement would 

i 

lower dust temperatures by a factor A 4+ £ (where N is the factor 
of emissivity enhancement), since the dust grains would radiate 
more efficiently. Coagulation could therefore explain both the 
mass discrepancy and the low dust temperatures we observed 
when fitting the spectral energy distributions with an emissivity 
index f3 — 2. 

IVoshchinnikov et al] d2006l) analyzed the influence of grain 
porosity on the dust properties. They show that as the poros- 
ity increases, the millimeter mass absorption increases and the 
dust temperature decreases. Using their results, we see that the 
mass absorption coefficient for a porous composite grain (amor- 
phous carbon+astr onomical silicate) can be up to 22 times the 
ILi & Draind J2001) value at 1 mm. Therefore, in theory the mass 
discrepancy between the virial theorem and dust estimates could 



be due to a larger fraction of porous grains or a larger porosity 
of the grains in the SMC GMCs than in the Milky Way. 

However, since grains coagulation is proportional to the den- 
sity of dust particles, one would expect less coagulation of dust 
grains in the SMC than in our galaxy, due to the lower dust con- 
tent. Indeed, recent models of the regulat ion of the grain size dis- 
tribution by shattering and coagulation (Hirashi ta & Yanll2 008) 
sugge st that, for a MRN grain size distribution dMafhis et al.l 
119771) . the importance of both these processes is significantly 
reduced at low metallicities like the one in the SMC. But the 
grain size distribution in the SMC seems significantly different 
than the one in our galaxy: the differences between the SMC and 
Milky Way diffuse ISM extinction curve and infrared spectral 
energy distribution indicate a larger fraction of the dust mass in 
small dust grains JWeingartner & Draindl200"lt iBot et al J 120041: 
iBernard et al.l [2008F 7 This enhancement in small dust particles 
is therefore likely to also impact the grain coagulation process. 
Using the dust grain size distributi on for the SMC and the Milky 
Way (for R v = 3.1) given by IWeingartner & Draind (2001) 
and a simple coag ulation model (using the collision rates due 
to turbulence from lDraind ( 1 1985b and assuming that grain stick 
when they collide), we observe that for a given total dust mass, 
small grains (10A < a < O.l/i m) could coagulate twice more 
quickly on big grains than in our galaxy. The enhancement of 
small grains in the grain size distribution of the SMC therefore 
seems to favor coagulation processes while the low metallicity 
of the SMC would penalize it. More detailed models would be 
needed to quantify the probability of large porous aggregates in 
the SMC. Such models are beyond the scope of the current paper. 
It seems unlikely at present that large porous aggregates would 
be much easier to form in the SMC than in our galaxy. But the 
type of aggregates formed in the SMC (and hence their porosity) 
could be different than the one in the Milky Way due to the dif- 
ferent size distribution and this could also potentially affect the 
dust emissivity. 

We can not completely exclude at present that the mass dif- 
ference we observe between the two methods (virial theorem 
and dust millimeter emission) is due to a higher sub-mm emis- 
sivity than that measured in Galactic GMCs. The presence of 
large porous grains in the SMC through coagulation would in- 
deed explain both the mass difference and the low temperatures 
observed. But our current understanding of dust in the SMC 
does not favor more coagulation in this environment than in our 
galaxy. 



5. Conclusion 

New LABOCA observations at 870yum were obtained in the 
south-west region of the Small Magellanic Cloud. All the GMCs 
observed previously in CO are detected. We build spectral en- 
ergy distributions by comparing the fluxes measured at 160/im, 
870//m, 1.2 mm, 3 and 6 cm for each cloud. We observed that 
the sub-millimeter and millimeter emission, in all clouds of this 
study, is always dominated by dust emission. Fitting the spectral 
energy distribution with a modified blackbody with a standard 
/3 — 2 gives surprisingly low dust temperatures Tdust ~ 12K. 
Given the uncertainties, we also took f3 = 1 (giving fitted tem- 
peratures around Tj usl ~ 19K) to deduce lower limits on the 
molecular cloud masses obtained from the dust emission. Using 
a dust emissivity at 870/mi that we determine in molecular re- 
gions of our galaxy and scaled to the dust-to-gas ratio of the 
SMC, mass estimates for the molecular clouds observed (as well 
as lower limits on these masses) are derived from the 870jum 
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dust emission and compared to virial masses deduced from CO 
observations. 

We confirm the discrepancy between the virial masses and 
the one deduced from dust emission that was observed previ- 
ously in the molecular clouds of the SMC. This discrepancy has 
now been observed at two different wavelength with two differ- 
ent instruments and as such can no longer be attributed to instru- 
mental effects. Furthermore, uncertainties on the origin of the 
sub-millimeter emission have been assessed (free-free contribu- 
tion, CO line contribution). The mass ratio M V j r /M"™ observed 
in the SMC is clearly below one, even for conservative assump- 
tions on the dust emissivity, dust temperature, etc. In contrast, 
for similar clouds in the solar neighbourhoods the observed val- 
ues are systematically above 1 . The direct interpretation of this 
result is that the CO line widths do not trace the full velocity 
dispersion of the gas in giant molecular clouds in the SMC. This 
could be due to the turbulent nature of the intersellar gas and the 
formation of CO in the densest structures. 

This low mass ratio in the SMC is unexpected for large 
molecular clouds and will have to be studied further in order to 
assess its origin. The dust properties in the SMC could be differ- 
ent than in our galaxy and we can not completely exclude that an 
enhanced grain coagulation in the SMC could increase the dust 
emissivity and lead to an overestimate of the true mass by the 
dust millimeter emission, Herschel observations will give more 
informations in that respect by complementing the wavelength 
coverage. The cloud+envelope geometry could also play a role 
in enhancing or decreasing the mass discrepancy and higher res- 
olution observations will be needed to assess this uncertainty. 

Further investigations are necessary in order to better under- 
stand the origin of this discrepancy, whether it is a general prop- 
erty of clouds in low metallicity galaxies and what it tells us 
about the molecular clouds properties. This will be better done 
with Herschel, which will be sensitive to dust emission on all 
angular scales and will provide a better spectral sampling of the 
far-infrared and sub-mm emission. Also, spectroscopic observa- 
tions of the C + line with Herschel will probe the extended en- 
veloppes of the molecular clouds and enable us to test whether 
the C + linewidth are comparable to the one observed in CO (i.e. 
test whether CO velocities trace the motions inside the molecular 
clouds). 
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Appendix A: Systematics 

A. 1. Effect of the local background subtraction on the 
masses derived from the sub-mm emission 

For this study, we tried to minimize the effect of CO-devoid en- 
veloppe on the dust emission by performing a local background 
subtraction. This option was chosen to be conservative in or- 
der to compute masses from dust sub-millimeter emission that 
are comparable to virial masses obtained from CO observations. 
However, it can lead to complicated uncertainties as discussed in 
Section |4~T1 

In this appendix, we test the impact of the background sub- 
traction on our results by redoing the analysis on the molecular 
clouds in the south-west region of the SMC, but without local 
background subtraction. Instead, we took the same background 
region for all clouds in all maps at the location 00:46:02.85 - 
73:19:26.5 (J2000). The fluxes obtained from 160yum to 6 cm 
are given in Table I A. fl On average, they are 1.5 times higher at 
870//m. All further steps of the analysis were then reproduced 
in the same manner as described in Sect. 12. 41 and the masses ob- 
tained are given in Tab. IA.2I Doing so, we find a median mass 
ratio M vir /M'™" = 0.14 (0.23 for the upper limit on the dust tem- 
perature), i.e. 1.5 times lower than with the local background 
subtraction. This confirms that our attempt to minimize the effect 
of an extended enveloppe by subtracting a local background to 
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Table A.l. Same as Tab.[T]but the radii correspond to deconvolved CO radii and the fluxes have been extracted without the local 
background subtraction. When the observed radius of the cloud was lower than the beam, the deconvolved radius is then set to NaN 
and we use the observed CO radius to compute masses instead. . 



name 


Radius 


S 160fim 


S ilQfim 


S l. 2mm 


S 3cm 


S 6cm 






pc 


Jy 


mJy 




mJy 




mJy 


mJy 




LIRS49 


10.6 


39.5 ± 3.4 


2241 + 266 


457 


+ 91 


36 ±5 


36 ± 


4 


LIRS36 


13.4 


24.5 ± 4.6 


1551 ± 140 


462 


+ 90 


80 ±30 


74 ± 


21 


SMCB1#1 


4.8 


3.9 ±0.9 


507 ± 


120 


215 


+ 63 


< 10 


< 5 




SMCB1#2 


9.5 


10.1 +2.1 


1252 + 185 


265 


±71 


21 ±6 


24 ± 


3 


SMCB1#3 


NaN 


5.4+1.4 


589 + 


99 


69 + 


19 


12 ±4 


13 ± 


3 


Hodge 15 


12.5 


0.68 ± 0.35 


467 ± 


83 


176 


±43 


2.3 ± 1.2 


5 ±2 


SMCB2 South 1 


8.4, NaN, 7.2 


45 + 9 


2603 + 557 


917 


± 187 


114 ±20 


105 ± 25 


SMCB2 North 2 


14.8, NaN 


43 + 6 


2137 : 


t330 


577 


± 165 


79 ± 15 


87 ± 


10 


SMCB2#5 


7.6 


6.3 + 1.5 


170 ± 


116 


71 + 


44 


27 ± 17 


27 ± 


17 



the emission observed in each cloud, decreases the mass discrep- 
ancy observed the submillimeter emission and the virial masses. 

A.2. deconvolved CO cloud radii 

We have used the sizes of the CO clouds as reported by Rubio et 
al. 1993b and Lequeux et al. 1994. Their sizes are the observed 
sizes which were not de-convolved by the beam size as the CO 
clouds most often showed larger sizes than the SEST beam. If 
one de-convolves their observed sizes assuming the standard 
aproximation of a source size smaller than the beam size, then 
the cloud sizes become smaller than the observed sizes published 
and used in the main text of this paper. The raddi obtained for the 
deconvolved CO clouds are given in Table lA.fl For a 43"arcsec- 
ond beam the de-convolved sizes are smaller than the observed 
sizes by a fraction of 0.3 to 0.8. In some cases, when the CO map 
had been fully sampled, the observed CO radius can be already 
smaller than the beam size. In this case, the deconvolved radius 
is marked in Table lA.fl as 'NaN' and we use the observed radius 
to compute virial masses. Since deconvolved radii are smaller 
than the observed ones, this results in smaller virial masses for 
the clouds by a similar fraction (see Table IA.21 . Thus the dif- 
ference between the gas mass derived from the dust emission 
becomes even larger. 

Combining the effect of deconvolved radii with the lack of 
the enveloppe subtraction (from Appendix lA. U . we obtain a me- 
dian mass ratio M,,,>/M™ m = 0. 10 (0. 16 for the upper limit on the 
dust temperature), i.e. twice lower than the one presented in the 
main text of this paper. This confirms that our different attempts 
to minimize known caveats like the presence of a CO-devoid en- 
veloppe and the resolution of the CO observations, decrease the 
mass discrepancy. The mass discrepancy observed in the south- 
west region molecular clouds could then be at least twice larger. 

Appendix B: Sub-mm dust emission outside the CO 
clouds 

The LABOCA 870pm map of the south-west region presented 
in this study shows extended sub-millimeter emission well out- 
side the molecular clouds observed in CO. In this appendix, we 
present a short analysis of this extended dust emission outside 
the CO peaks. We build the spectral energy distribution from 
160pm to 6cm by extracting fluxes at each wavelength for a 
region defined as being observed at all wavelengths and be- 
ing outside the CO peaks studied. This region hence defined 
corresponds to a solid angle of 1.00787 x 10~ 7 sr. For the ex- 
tended emission observed in this region, the spatial filtering in 
the SIMBA 1 .2mm observations is no longer negligible and the 



Table A.2. Gas masses deduced from the sub-millimeter dust 
emission (as Tab. [3] but without local background subtraction) 
and virial masses deduced with the deconvolved CO radii. 



name 


H 


(12 K) 


1Y1 H 


(19 K) 


M yir 




10 4 


M 


10 4 


M Q 


10 4 MO 


LIRS49 


75: 


±9 


50 


±6 


7.7 ± 


0.5 


LIRS36 


57: 


±5 


34 


±3 


3.7 ± 


0.5 


SMCB1#1 


26: 


±6 


12 


± 3 


0.9 ± 


0.2 


SMCB1#2 


56: 


±9 


28 


±4 


4.3 ± 


1.4 


SMCB1#3 


21 : 


±4 


13 


±2 


2.5 ± 


0.2 


Hodge 15 


39: 


±7 


11 


±2 


5.0 ± 


0.4 


SMCB2 South 


95 i 


:21 


57 = 


t 13 


9.5 ± 


2.3 


SMCB2 North 


72 i 


: 12 


47 


± 8 


7.8 ± 


1.9 


SMCB2#5 


4h 


:3 


3 = 


t 2 


2.6 ± 


0.3 



flux obtained at 1 .2mm is considered as a lower limit. The spec- 
tral energy distribution obtained is shown in Fig. IB. II 

A free-free emission model is fitted to the 3cm radio flux 
and extrapolated to 160 and 870pm. The free-free contribution 
is then removed from the observed flux at 160 and 870pm be- 
fore fitting a modified black body spectrum with a spectral index 
of 1 (blue dashed line) or 2 (plain red line). The dust tempera- 
tures deduced are 

T dult = 14.8Kand 7^ = 10.8K, respectively. 
The dust temperature observed in the region surrounding the CO 
peaks is therefore similar to those observed in the CO-detected 
regions in our observations. 

The hydrogen mass of this region deduced from the dust 
emission fit is 9.3 x 10 4 M© for a standard spectral index /3 = 2 
(5.2xl0 4 MO for/3 = 1). This is comparable to the mass, deduced 
from the sub-mm dust emission, of the smallest CO-detected 
molecular clouds of this study. 

Using the Hi observations of the SMC from 
IStanimirovic et all (1 19991) . we estimated the neutral hydro- 
gen mass in this region to be 2.5 x 10 4 MO, i.e. 2 to 4 times lower 
than the mass of the gas deduced from the dust emission. Most 
of the gas observed through the extended sub-mm emission 
outside CO-detected molecular clouds is then probably cold 
molecular hydrogen. This is consistent with the fact that we 
observe no correlation between the 870pm emission and the Hi 
21 cm emission. 
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Fig. B.l. Spectral energy distribution from the far-infrared to the 
radio, observed for the region outside the molecular clouds ob- 
served in CO. The fluxes observed at 160, 870/im are fitted with 
a single modified black body either with an emissivity index of 1 
(dashed blue line) or 2 (plain red line). Radio fluxes computed in 
the same regions are overplotted and the extrapolated free-free 
emission is displayed as a black line. 



